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1 Introduction 



Over the recent years, the gauge-string duahty has proved to be a valuable tool for under- 
standing strongly coupled conformal field theories. In particular, the shear viscosity t] of 
a wide class of strongly-coupled gauge-theory plasmas has been computed. These include 
the A/" = 4 super- Yang-Mills theory [H [21 [3], the theory dual to the near horizon limit 
of M-branes [H [5], some non conformal theories [6l [71 [H [9l [10], non relativistic theories 
[T2l [T3] and extensions of these theories which include the introduction of a chemical 
potential [TU [13 [HI [IT], fundamental matter [TSl [TH] and non-commutativity [20]. In all 
the cases considered the large 't Hooft, large value of the ratio of the shear viscosity to 
entropy density was 

s 47r ^ ^ 

General arguments showing that the ratio ([H) is universal can be found in [21], [22l [23] . 
Extensions of (HD to finite 't Hooft coupling or can be found in[ai[25l[26l[27l[2Hl[29l[3ni[3T]. 
Perhaps of more practical importance is the proximity of ([1]) to the apparent shear viscosity 
to entropy ratio of the quark gluon plasma, presumably observed in heavy ion collisions 
[32| [33| [M] . Discussions of the universal behavior of other properties of the quark gluon 
plasma can be found in [33 [Ml [33 [3B] • 

To leading order in the hydrodynamic approximation, in the absence of generalized forces 
and when there are no spontaneously broken symmetries, the shear viscosity is the single 
parameter which is required in order to completely specify the stress-energy tensor in a 
conformal theory on R'^'^E Indeed, when the mean free path of the theory is smaller than 
the typical momentum scale, one can parameterize all the quantities appearing in the hy- 
drodynamic equations by the hydrodynamic variables u'^, e, and fia- The velocity field m'^ 
is defined as the velocity needed to boost a fiuid element to its rest frame, e is the local 
energy density in the rest frame of a fiuid element, and the fia are local chemical potentials 
conjugate to charge densities pa also defined relative to the rest frame of a fiuid element. 
With this parameterization, conformal invariance dictates that 

{Tf,u) = ■j^eidUf.u^ + 7]^y) - 7]a^y + 0{d'^) , (2) 

^ As summarized in |39| , a broken symmetry implies an extra velocity field associated with the coherent 
motion of the condensate. With the appearance of an extra velocity field, one should be able to construct an 
extra viscous term in the energy momentum tensor. A holographic construction of a two-component fluid 
can be found in [30] . 
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where 

a^,y = 2d^^Uy> (3) 
and triangular brackets denote a symmetrized, traceless projection: 

A<,.> = \P,^P/ (Aa. + A^a) - ^P^.P'Ma. (4) 

with 

P^iu = V^lu + u^Uy (5) 

a projection onto the space orthogonal to the velocity field. 

The 0{d'^) in (I2j) refers to expressions which are subleading in gradients of the hydrody- 
namic fields. I.e., they include at least two derivatives of m^, e and ^a- At second order in the 
hydrodynamic expansion many more parameters are needed to completely specify the stress 
energy tensor. We call these parameters second order transport coefficients. Four of these, 
and Aj with i = 1, . . . , 3 were discussed in [H]. They are associated with the contributions 

of 

S^l^ = a<,.A^\) , = a<^A^\) , Sjfj = u^^^Lo\ , (6) 

with 

uo,, = (9am, - d^ux) , (7) 

to the second order stress-energy tensor: 

{T,.) = ^e{du^Uy + r/^.) - r^a,, + r^r/S^ + + 0{d'). (8) 

The coefficients r^r and Aj were computed for theories whose dual is given by Einstein gravity 
[m SSI HSl HH HSl HH], five dimensional Einstein-Maxwell gravity with a single f/(l) field 
[171 HH] or five dimensional Einstein gravity in the presence of a dilaton source |19]. In ^7] 
it was observed that in all these examples, the transport coefficients r^, Ai and A2 satisfy 

4Ai + A2 = 2r]T^ (9) 

when the boundary theory spacetime dimension is larger than three. In three dimensions 
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S^i^ = 0, so it makes no sense to define the associated transport coefficient Ai. From now 
on we will always be considering d > 3. In this work, we show that iQ holds for a class 
of strongly coupled, conformal, gauge theory plasmas with conserved global charges which 
have a bulk dual. We assume that the bulk dual of the stationary fluid can be described by 
a stationary and static black hole and that the black hole solution can be extended in the 
sense of jl3]. Further, we show that A2 is given by 



^2 = / ^ dx , (10) 

J r-i- 



where s{r) and /(r) are given by the metric components of the bulk dual to a static, sta- 
tionary configuration of the fluid. 



d-l 



ds^ = -f{r)dt^ + ^dxl\+ 2s{r)drdt, (11) 



r_|_ is the largest zero of f{r), and is the {d + l)-dimensional gravitational constant of 
the bulk theory (see below). Equation ( ITOi) coincides with what was found in [16] where the 
case /(r) = 1 — (r^/rY, s(r) = 1 was studied. Note that (Q imposes a constraint amongst 
the transport coefficients of the theory, whereas ([1]) relates a transport coefficient (the shear 
viscosity) to a thermodynamic variable (the entropy density). In that sense ([T]) and ([9]) differ. 
Further discussion of and fllOj) and their possible extensions can be found in section [31 
The derivation of these two equations can be found in section [2l 

2 Computing (T^j,} 

Consider an action of gravity coupled to matter, 

1 f ~ .D-i 



S = -TTTTi I V^iR + ^m) d^'-'xdr + Scs- (12) 



2k 

The matter Lagrangian may include scalar s "J)"", and Abelian vectors Al^, 

Cm = -G^pi^)dM^''d''^P + \/($) - Kabi^) (13) 

with -FXfjv = "^M^AT — dNA.1,j. The indices a, b and a, P run over the number of gauge fields 
and the number of scalar fields, respectively. The indices M = 0, . . . , D — 2, D run over the 
D AdS directions with the D coordinate denoting the radial AdS direction r. The indices 
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yU,, V which have been introduced in the previous section, run over the d = D — 1 directions 
transverse to the radial AdS direction. These will be called the transverse directions for 
short. The action ^cs denotes a possible Chern-Simons term. The Einstein equations which 
follow from the action f[T2l) are 

Rmn = ^MN + Qmn^ , (14) 

where 

Tmtv = Gaf^dM^'^dN^^ + 2KabF^^F'j,'' , (15a) 
1 = --(D^) (^("^^ + K^b{^)F^MNF'''"') . (15b) 

In f fT^ the matter stress-energy tensor has been decomposed into a term Qmn^i whose index 
structure is proportional to Qmn^ and another term ^mn whose index structure comes from 
the gauge fields, their derivatives, and derivatives of the scalar fields. 

Since we are interested in a thermal state in the dual theory, we consider black hole 
solutions to ([HD. Let us assume that the equations of motion allow for a translation and 
rotation invariant (in the spatial transverse directions) time- independent, asymptotically 
AdS black brane solution. In the Eddington-Finkelstein coordinate system, the most general 
such solution takes the form fllip . The asymptotically AdS boundary is located at r ^ oo 
where f{oo) = 1 and s(oo) = 1. We denote the largest zero of /(r) by r+. This is the location 
of the outermost horizon of the black hole. Assuming that the gauge and scalar fields are also 
time-independent and invariant under translations and rotations (i.e., = $"(r). A"} = 0, 
A'^ = A'^{r)), we find that 

%i, = , Toi = , (16) 

where the indices i and j run over the transverse, spatial directions. To make contact with 
the notation in the rest of the paper we denote the components of %mn and T by T^/*^, 

Inserting fITB]) into f[T^ we can determine /(r) and s{r) in terms of X'-"^ and X^]^ through 
the M = i, N = i and M = D, N = D Einstein equations. We find that 

s(r) = e-^^'°°^^ol,(^)««^ (17b) 
One can show that these solutions are consistent with the remaining Einstein equations as 
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long as the bulk energy momentum tensor is conserved, D^T = 0. The extra integration 
constant, (e — e^) /(2k,'^(D — 2)), has been chosen to be the coefficient of the r~^^~^^ term 
in a series expansion of /(r) near the asymptotically AdS boundary. The reason for this 
notation will become clear shortly. The requirement that the boundary is asymptotically 
AdS has fixed the upper end of the integral in fll7bl) . ensuring that s{oo) = 1. It also implies 
that %^^\oo) = —{D — 1) in order that /(oo) = 1 be satisfied. 
In this work we will further require that 

/(r) - 1 = C(r-(^-^)) , s{r) - 1 = C(r-(^"^)). (18a) 

The reason is as follows. The bulk equations of motion specify the metric up to two sets of 
integration constants. The ffist set specifies the value of the metric at the asymptotically 
AdS boundary which in turn specifies the metric of the dual CFT. The second set of integra- 
tion constants may be associated with the 0{r~^^~^'>) coefficients of a near boundary series 
expansion of the bulk-metric (as specified above) and are associated with the energy mo- 
mentum tensor of the dual CFT. Once these two sets of integration constants are specified, 
the full bulk metric is completely fixed. Its near boundary Taylor series can be determined 
by iteratively solving the equations of motion. Generically, for odd D a non-fiat boundary 
metric induces logarithmic terms in the near boundary series expansion of the bulk metric. 
At order (9(r~*^^~^Mnr) these logarithms are associated with conformal anomalies which 
render the trace of the boundary theory stress-energy tensor non-vanishing (as expected 
when the boundary theory geometry is non-flat). But it is not only the boundary metric 
which can induce such logarithmic terms. In odd dimensions, any time the bulk stress energy 
tensor generates an (9(r~^™) term in the metric, with m < {D — l)/2, there will also be an 
(9(^7--(^-i) inr) term which implies that the state in the boundary theory is not conformally 
invariant. An explicit example of such behavior can be found in [SU]. We will also restrict 
ourselves to solutions for which 

a; = o(r-(^-^))5; (18b) 

and will set the sources dual to the scalars to zero. An example of an asymptotically AdS 
background satisfying flTSl) is given in [211 We discuss the effects of loosening the 

restriction (fTSi) in section [3l 

^Note that we do not require the gauge field to vanish at the future horizon and therefore it is hkely that 
it diverges at the past horizon. In fact, it is likely that the whole perturbative solution diverges at the past 
horizon because generic solutions of viscous fluid dynamics are not expected to be regular in the infinite 
past. We thank A. Karch, D. Son, and especially R. Loganayagam for clarifying this point. 
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We may now use the prescription of [531 EH ESI ES], adapted to Eddington-Finkelstein 
coordinates to convert the metric (ITT!) to the boundary theory stress-tensor (T^u)- We 
find that 

{T,u) = ^ {dSy. + V,u) + ^ {dSX + V,u) . (19) 

The first term in ( |T9l) represents the component of the energy momentum tensor associated 
with the extrinsic curvature of the near boundary metric of the black hole and a boundary 
cosmological constant counterterm to the bulk action. The second term in f|T9|) represents 
extra contributions arising from the matter content of the theory, and from possible matter 
counterterms when holographically renormalizing the bulk action [ST]. Since we have as- 
sumed that the theory is exactly conformal, these extra contributions can be determined up 
to an overall factor, eo. 

So far, we have considered a stationary and static black hole. Following the philosophy 
of [13], we assume that the black hole solution (ITTl) can be extended to capture the fiuid 
dynamical nature of its boundary dual. We start by boosting the transverse coordinates by 
a constant velocity u'^. After this boost, the metric ffTTj) takes the form 

ds'^ = {-f{r)Uf,Uudx^dx'' + P^^dx^'dx") - 2s{r)u^dx^dr. (20) 

Clearly, fl20!) is still a solution to the Einstein equations since it has been obtained from 
the solution ffTTj) by a coordinate transformation. The stress energy tensor of the boundary 
theory dual to (1201) is the boosted version of the stress energy tensor f[TI?l) 

{Tf,u) = {dUf,u^ + Vf^u)- (21) 

As explained earlier, in the boundary theory, the hydrodynamic approximation involves 
allowing the velocity field, energy density and chemical potentials to vary slowly in the 
transverse directions. Thus, to obtain the gravity dual of hydrodynamics, we allow the bulk 
quantities, m'^, e and the various charges to slowly vary with x^. Once we do that, (l20l) 
will no longer be a solution to the equations of motion. Therefore, we add a correction g^^J 
to the metric (1201) . linear in derivatives of the velocity fields, energy density, and whatever 
chemical potentials are present. Similar corrections, and are added to the scalars 

and gauge fields. The explicit values of gj^) , and A'^^^^ are determined by solving the 

Einstein equations (I14p together with the equations of motion for the matter fields. Once the 
corrections which are linear in the derivatives are obtained, one can continue this procedure 
to obtain two-derivative corrections gjul, and then three-derivative corrections. 
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etc. Such a perturbation scheme has been shown to work (up to the two derivative level) 
in the case of Einstein gravity in various dimensions [HI 1131 SH SSI SS] , in the case of 
Einstein theory coupled to a dilaton [19] and Einstein-Maxwell gravity in five dimensions 
with and without a Chern-Simons term [13 HHl [58] . 

After carrying out this perturbative scheme, one will end up with a (perturbative) solution 
to the Einstein equations dual to a non trivial fluid flow in the boundary theory. In [IH] 
it has been shown that conformal invariance of the boundary theory implies that the bulk 
metric must take the form 

ds^ = -2S (^dr + r (^u^dxu^ - 'J~^^^^^^''^ ^^"^ Uf.dx^' 

+ r'^Ku^uJx^dx" + r^P^^dx>'dx'' + r^Jx {P\u^ + P\u^) dx^dx'' + r'^Tr^^dx^'dx" , (22) 

where S{r), K{r), J^{r) and vr^i/(r) are invariant under the boundary theory Weyl-rescaling 
^7^1/ G~'^^ri^y and depend implicitly on x^jfl In addition, the traceless tensor is orthog- 
onal to the velocity field u^. If we set 

5(r) = s(r), ir(r) = -/(r), J^(r) = , 7r^.(r) = , (23) 

and set to a constant, we recover our boosted stationary solution fl2Ul) . Information on 
the viscous corrections to the boundary theory energy momentum tensor is contained in n^^y. 
This follows from the mapping of the bulk metric to the boundary stress energy tensor. In 
our notation (see 



{Tf,u) = -r r (dUf^u^ + rjf,^) —Finite lim r drTT^^{r) + . . . , (24) 

a — i Z Lr^oo J 

where Finite [limr^oo -B(^)] means the finite part of B{r) as r is taken to infinity. The ellipses 
in (12^ represent corrections to the Brown York stress tensor coming from holographic renor- 
malization and from the matter content of the theory. Thus, to find the various transport 
coefficients associated with the boundary theory fluid, we need to compute vr^^y and deal with 
the various corrections to (T^;/) coming from holographic renormalization and from varying 
the matter action with respect to the boundary metric. 

As a warmup, let us compute the shear viscosity of gauge theory plasmas whose dual bulk 



In [49] the bulk form of the metric was shghtly more general than l|22|) in that the coefficient of the 
r^Piiu term was an additional Weyl invariant scalar. Here we have chosen a gauge where this extra scalar 
has been set to one. 
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action is of the form (fT2!) . For this we need to find the stress energy tensor of the plasma 
to first order in the hydrodynamic approximation. Equivalently, in view of (1241) . we need to 
find vr^i/ to first order. As stated earher, we denote this term and in what follows we 
will use the notation 5^"^ to specify the n'th order component in a derivative expansion of 
the field B. 

As pointed out in [13], a derivative expansion of the Einstein equations results in a 
differential equation in the radial coordinate only. This implies that it is sufficient to solve the 
Einstein equations in the neighborhood of an arbitrary point Xq in the transverse directions 
for all values of r and then extend the solution to the entire manifold using Lorentz invariance. 
A further simplification can be made if we choose a coordinate system where 



m''(xo) = (1,0,...,0). (25) 

This is a simplification since in the neighborhood of Xq, the Einstein equations naturally 
decompose into tensor, vector and scalar modes of the SO{D — 2) symmetry which leaves 
u^{xo) invariant. The equations of motion for the tensor modes, 7rjj(xo), can be read off from 
the M = i, N = j^i components of the Einstein equations (|T^ expanded around Xq, 



I^ = c^^+^^+^. (26) 

Working out R^^ explicitly, using gf^ = X-j^ = for i ^ j and using (I17ap we find that fl26l) 
reads 

'^vr;fy=-(Z^-2)r^-V.,, (27) 

where a prime denotes a derivative in the radial direction. After extending the solution to 
f l27|) to the entire manifold it takes the form 

7rW(r)=a^,F(r) (28) 

with 

r {^""^^ - rf-')s(x) , , , 

While vr^i/(r) can only be written down implicitly, we can use ( !29l) to obtain 



Finite 



lim r^9,7r(^j(r)l =r:;^-2. (30) 



In order to use to evaluate the boundary theory energy momentum tensor, we need 
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access to the extra terms in arising from holographic renormahzation and the matter 
fields. By using fll2p and fllSI) and a simple power counting argument, one can show that the 
gauge fields do not contribute to such terms. We elaborate on this argument in appendix 
1X1 We conclude that only terms involving derivatives of the scalars can contribute to the 
boundary theory energy momentum tensor. As discussed in the introduction, for the theories 
we have in mind, such expressions can not contribute to {T^u) at the one derivative level. 
Thus, 



To compare this with ([1]), one can compute the entropy density of the boundary theory fiuid 
from the Bekenstein-Hawking formula: 



47rr^ 

■ 



(32) 



The ratio rj/s agrees with ([T]). For the class of theories we are considering, this result was 
previously obtained in |26j . 

(2) 

We proceed to compute ttIiJ , the second order contribution to tt^j^. The computation is 
similar to the one that lead to fl2Sl) . At second order in the derivative expansion, ttIiJ can be 
computed from the i j components of the Einstein equation expanded around Xq 

1^ = ^^ + ^^+^. (33) 
Writing this explicitly in terms of the metric components, we find 



^ ^ n=0,...,3 71=1,2 



(2) 



o t(1)\ om 0^^£'-3 t(1)o „, D'^ j) , ^ ^ U 



- 2 d^J^ - 2{D - 2)r^-V«9o«^.^ - r 



s s 
- 2r^-%d,ys + ''''~'^(^'^^)' + 2r (r^"^9(,.9o«,))' . (34) 

The quantities P„, P„ and T are functions of r, s(r), /(r), S^^\r), K^^\r), F(r) and their 

— (n) 

derivatives. We refrain from writing them out explicitly. The S^^, 's are given by 

^t^j = d^ui^^d'^Uu) , S^^^j = u^di^^d^Uu). (35) 
The superscripts of the S's and S's are a means to enumerate these components and are 
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not related to the hydro dynamic expansion. is the first order correction to defined 
in fl22l) . expanded around Xq. It can be determined in terms of the vector components of 
'^MTV t>y solving the first order vector mode equation equivalent of ( l26l) . The expressions 
for S^^^ and K^^^ can also be determined by solving the first order scalar equations. It is 
straightforward to show that S^^^ = 0. In (IMI) terms proportional to have been omitted. 

Notice that the index structure of vr)- is inherited from the index structure of the terms on 
the right hand side of (1341) . Further, the index structure of the boundary theory stress energy 
tensor is equal to that of VTjj plus any terms coming from the matter fields or counterterms 
from holographic renormalization. Before proceeding, let us take a moment to characterize 

(2) 

the possible terms which can contribute to vr^^ . Recall that ttij is a Weyl invariant traceless 
symmetric tensor orthogonal to the velocity field. At second order in the derivative expansion 
there are a finite number of combinations of derivatives of e, and which can contribute 
to irlj . If we consider terms composed only of derivatives of and e and use energy- 
momentum conservation, d^iT^J) = 0, or 

P/d^e = -edu^'d^u^ + 0(9^) , (36a) 
u^d^e = -e-^d.u" + 0{d'), (36b) 

then the only Weyl invariant, even-parity terms which we can construct are H^^J with 
n = 0, . . . , 3. If we loosen the restriction on Weyl invariance there are two more possi- 
ble terms: S^^^j and E^f] [12] • Allowing for pseudo-tensors, or terms which involve /i^ and 
their derivatives, will lead to many more possibilities. We denote these other terms by 
H^^J with n > 4. We exclude from this list expressions which are equivalent under current 
conservation = 0, or 

u^'d^Pa = -Pad.u'' + 0{d^). (36c) 

In what follows, we will rewrite all the terms on the right hand side of (l34l) in terms of the 

basis expanded around Xq. Our assumption regarding conformal invariance dictates 

that the coefficients of the non Weyl-invariant expressions will eventually vanish. However, 

— (") 

as we will see below, it will be useful to keep track of the coefficients of the S^^, terms. 

Let us start with the (T — 2r^~^sFT^^'') aij term in the first line of ( !34|) . Using ( !36l) . at 
the one derivative level we find that the only scalar which can be constructed out of m^, /Zq 
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and e is d^u'^. Thus, 

(r(r) - 2r^~hF%^^\r)) = e^{r)d^u^'. (37) 

Next consider the various derivatives of J^^^ . The Weyl invariant vector or pseudo- vector 
jIP is constructed from m'^, e and /Iq and one derivative. From (!22l) we see that only the 
component of J^^"* which is orthogonal to the velocity field will enter the metric. Thus, J^^ 
must be composed of 

a.^, d,^ (38a) 
e //ft 

and, in d = 4, 

= e^.x.u^'d^u^ (38b) 
Using fl36al) . we find that the expressions involving J^^^ in 0341) must take the form 

02{r)doU(^jdoUi) + B^{r)di^i {u"d^Uj)) + . . . , (39) 

where the ellipses denote expressions which can not possibly contribute to E^-"^ with n <2. 
(Note that £^ is a pseudo-vector and that ^{^^v) = 4S|fj].) Like \ s is also Weyl invariant. 
Therefore, the third line in flMl) will also take the form fl39l) . 

We are left with the T\j term in the first line of flM|) . defined in fll5ap . This term has 
two components, one involving the gauge fields, and the other involving the scalars. Using 
the same argument as in the last paragraph, we find that the scalars can contribute to the 
02 term in (!39|) . To see what the contribution of the F^j^F\^ term in fllSaP looks like, we 
expand the gauge field A'^ in a derivative expansion, i.e. 

= c'^(°)n^ + (P;<(^) + c«(^)n^) + . . . . (40) 

The first term on the right hand side of (HOl) . c°'^^\ is the zero order stationary solution 
boosted to a velocity u^. The second expression, ^P^a"''^^ + c°'^^^Uf^, is the first order cor- 
rection to separated into a contribution orthogonal to the velocity field and a contribution 
parallel to it. As was the case for {T^^), in order for the boundary current (J^) to trans- 
form homogeneously under Weyl rescalings, we need that at^^^ be a Weyl invariant vector. 
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Plugging (HOj) into F^j^F\^ and expanding around Xq, we find that for i j 



.4^-2^a(0)^6(0)v(3) + !^/«;'a)'«Ml)' 



-.2 « J 



a(l)/ 



- (41) 

We liave omitted terms wliicli depend explicitly on x'^. The first expression contributes only 
to S^^^ and the remaining terms can only contribute to 62 in fl39|) . or to H^f^J with n > 3. 
We can summarize the previous discussion by recasting flM|) into 

/ n=0,l,2 n=l,2 



(42) 



{n) 



where the ellipses represent expressions which do not contribute to S)^^'' and s|t"^ with 



n < 2. The last three terms in fH2]) are symmetric traceless tensors composed of the ve- 
locity fields and their derivatives. Therefore, they may be expanded in terms of the basis 
{Sj^°\ . . . , Tif^ , E'J'', ^fj}. Doing so, equation (142!) takes the form 



^ n=0,l,2 n=l,2 



(43) 



with 



Ao(r 
Ai(r 
A2(r 
Ai(r 
A2(r 



i'o(r) + 
Pi{r) - 

P2{r) + 
Pi(r) + 
P2{r) - 



3{D 


-2) 


D- 


f4 


3{D 


-2) 


2{D 


+ 4) 


6{D 


-2) 


D- 


f4 


6{D 


-2) 


D- 


f4 


6{D 


-2) 


D- 


f4 



^i(r 
^i(r 
^i(r 
^i(r 
^i(r 



+ 



D + 4 
D-2 
4(D + 4) 
D-2 



+ 4 
D-2 

D + 4 
6 

L) + 4^ 



02{r 
02{r 

e2{r 

doir 



Mr 



D + 4 
^-2 ^ , 



4(D + 4) 
6 



+ 



D + 4 
6 

D + A^ 
D-2 



(44a) 
(44b) 
(44c) 
(44d) 
(44e) 



_____ fn\ 

So far (H3l) shows us the contribution of vr^^- to the boundary theory stress tensor. Further 
contributions to {T^y) will come from the matter action and from holographic counterterms. 
We have argued previously that these can only come from the scalar sector of the theory. 
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The only possible contributions from this sector which are relevant to our analysis are trans- 
verse, traceless derivatives of the Weyl invariant quantity /i^/e. We can use the same logic 
that brought us from fl38|) to fl39l) to argue that these contributions will be proportional to 
dQUi^jd^Ui) and 9(j {u'^dfjUj^^. Now, since the 6i terms in (jUj) are (and will remain) unspeci- 
fied, we can work with an effective 6*2 and 9^ which capture the contributions of these extra 
terms from the scalar sector. In what follows we assume that such a procedure has been 
carried out, and then the sole contribution to r^, Ai and A2 comes from the Aj, i < 2 defined 
in (144|) with unspecified 6i. 

We are down to solving (143|) . Currently, this is made difficult due to the unknown 
functions 6'j. Fortunately, we can do without them. Recall that the coefficients of the S^^ 
terms must vanish if the boundary theory energy momentum tensor is to be conformally 
invariant. This implies that 

- + 6^^ + 6^:^^^^ ' (45) 

Substituting (H5!) into (H3l) and inserting the explicit values for the P„'s and P„'s, we find 
that 

Ao(r) = r (r-°-^s(r))' - 2r^ [r^F{r)^ + ^^3(r) , (46a) 
Ai(r) = Qr^-3s(r) - (r^-^ _ ^d-2~^ ^(^)^ ^ l^^(^) ^ (45^) 

A2(r) = -2r^~^s(r) - 4r V F{r)^ , (46c) 

and, of course, the A^'s vanish by construction. From fH6|l it follows that 

4Ai + A2 - 2Ao = 4 ((r^-2 - r^-^) F{r))' (47) 

and using flM|) and ([8]), we obtain the proposed relation ([9]). Since A2 in fl46cl) does not 
depend on ^3 we may use fl2^ and ([8]) to evaluate A2 explicitly. The result is given in ffTOj) . 
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3 Discussion 



In section ([2]) we derived our main result, which we reproduce here for convenience, 

4Ai + A2 = 2r]r^. (Q 

The transport coefficients A2, Ai and t^, are associated with H^^u and sjfj, respectively, 
cf. ([8]). Another result we obtained is an explicit expression for A2 in terms of the components 
of the bulk metric, dual to a static configuration of the fluid, cf. fllUp . A similar expression 
for Ai is 



Ai = rf "^s(r+) Finite 



4 



+00 



Sir 



lim I / 6^{x)dx 







(48) 



where, within the current analysis, 6*3 is an undetermined function 

The choice of basis {S^"''} for the second order contributions to {T^^) is arbitrary and 
was taken from [?T]. There, was chosen to conform with Israel-Stewart theory in the 
linear regime. But, one can, for example, work with the alternative basis 

y(0) _ y(0) , c.y(2) 

yW _ _ 4^(2) ^ ' 



which is similar to that of Using we find that in the alternative basis the coefficient 
of T^nu , A2, vanishes, 

r/r.sW + AiSW + AsS^^j = r/r.S(S + AiS«. (50) 

It is interesting to compare (Q with a weak coupling analysis. Such an analysis has been 
initiated in [H] and carried out explicitly for QCD (with and 3 flavors), QED and 0"^ theory 
in 4 spacetime dimensions and in the absence of a chemical potential in ^U\. According to 



^ The 03 term in originated in di^idj-^^ terms which showed up when we took derivatives of Weyl 
invariant quantities. It was the di^idj^e contribution to such second order derivatives which, when evaluated 
under the constraint equation (|36ap . generated a 6'3(r)9(j (u^S^Uj^) term. It seems hkely that by an appro- 
priate choice of the S]jjJ/'''s with n > 3, say, S^""'' = I?^^!?^^,^ with 2?^ the Weyl covariant derivative of [59] . 
one should be able to reshuffle the distribution of the O'i terms, and get rid of the 6*3 dependence in (|46p . 
The current formulation of the problem does not allow us to check this possibility explicitly. In the specific 
case considered in [JT], Ai was given by (|48|) with Oj, = 0. 

^ In [4^ the authors considered a CFT with a non flat boundary metric. This allows for an extra 
contribution to the energy momentum tensor at second order in a derivative expansion, not considered in 
this work, proportional to the Weyl tensor contracted with the velocity fields. The basis used in equation 
(1.3) of [46] and the discussion in [46] suggests an additional relation, similar to ([9]), for conformal theories 
on curved manifolds. 
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[m [60] the relation (Q does not hold at weak coupling. Instead one finds 



A2 = -2r/r,. (51) 

In the alternative basis, ([9j) and (15T|) can be translated into the statement that A2 increases 
from — 4(77X7^— Ai) at weak coupling to at strong coupling (and infinite N). The perturbative 
computation of [60] for QCD indicates that, in the alternative basis, A2 increases with the 
coupling, in agreement with our previous statement. 

To obtain iQ, we had to assume that the theory at hand is conformal, that its bulk 
dual can be truncated to an action of the form f[T^ . that the U{1) symmetries are not 
spontaneously broken, and that the various fields have an asymptotic expansion as in ([T8|) . 
The main reason why we worked with the asymptotic expansion in ( IT8|) is that, otherwise, it 
would be difficult to handle the extra terms which arise when holographically renormalizing 
the theory. 

Some of these restrictions could be easily removed. For example, we may allow for sources 
for the scalars which correspond to operators with dimension smaller than four. In such case, 
the main change in our computation is that the energy- momentum tensor may be sourced 
by gradients of the sources (as was the case for the dilaton in ^1]). This would simply mean 
that the right hand side of fl36a|) and (I36b|) needs to be modified by adding to it appropriate 



derivatives of the scalars. Since the merit of (I36ap and (136bp is in converting space-time 
derivatives of the energy density into spacetime derivatives of the velocity field, most of our 
analysis will remain unchanged. 

Perhaps a bigger challenge to ([9]) is a non-conformal theory. Recall that we relied heavily 
on conformal invariance to obtain the form of the first order vector expressions in f[38|) and 
to set the A„'s to zero in ( l43j) . If these coefficients were not zero, we would not have been 
able to obtain (H7I) . Curiously, even if the A„'s did not vanish, we could still have written 

4Ai + A2 - 2Ao - A2 = 4 _ ^d-2) ^^^^y ^ ^52) 

which suggests a relation of the form 

4Ai + A2 - 2r/r, - A2 = , (53) 
where the Aj are the second order transport coefficients associated with S^^, when expanding 
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a non conformal stress energy tensor in a hydro dynamic expansion, 

(T,,,) = eu^u, + pP^, - ria^, - CP^^.d^u^ + r/r.S ^ + A,S» + A.sJJ^ + . . . . (54) 

Of course, once the system is non-conformal there will be further corrections to coming 
from the bulk viscosity ( and pressure p. It remains to be seen if there is a relation of the 
form iQ or fl53|) that holds for non-conformal theories. Perhaps a good starting point for 
such an exploration are the relatively simple non-conformal backgrounds constructed in [61j . 
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A Holographic renormalization 

By varying the bulk action S'totai with respect to the boundary theory metric g^^J , one obtains 
the stress energy tensor of the boundary theory 

The action S'totai is composed of the bulk action ( fT2|) . evaluated on-shell, a Gibbons-Hawking 
boundary term, and a counterterm action 

Set = [ d'^xCct , (56) 

Jr=Ro 

which renders S'totai finite when evaluated on-shell. The contribution of the Einstein part of 
the bulk action and the Gibbons-Hawking term to (T^^) together with a boundary cosmo- 
logical constant was given explicitly in (!24|) . 
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Let us rewrite the counterterm action in the form 

Set = j d^x^iCo + C^^'A,^ + L^^^-A^^A^^ , + ...), (57) 

where the dependence of the Lagrangian on the gauge field has been made exphcit. The 
constructed from the various scalars and derivatives and can be thought of as 
acting on the terms on their right. We have denoted the metric on a shce of constant r = Rq 
by 7^jy. Near the boundary 7^^^ ~ i?o ^ ^7 ~ Rq- From 1^, A"^ R^"^^^ and assuming 
that we are not turning on relevant operators in the conformal boundary theory, 0" ~ R^°' 
with rria < 0. 

Suppose ~ R^. If the term is to have any contribution to the counterterm 

action in the i?o — ^ oo limit, we need that m — d{n — 1) > 0. On the other hand, £^i...;i„ 
is composed of the scalars and their derivatives. By our previous discussion this means that 
m < 0. Thus, only the Cq and £^ terms can contribute to the non- vanishing parts of in 
the Ro ^ oo limit. Gauge invariance and the absence of charged scalars imply that the 
term involves a derivative acting on A^. In this case involves three derivatives, so that 
C^Afj_ ~ Rq'^~^. Hence, this term also vanishes in the Rq ^ oo limit. We conclude that 
the contribution of the counterterm action to the boundary theory energy momentum tensor 
will not include the various gauge fields. 
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